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Previously developed "stochastic representation of deterministic interactions" enables exact treat- 
ment of an open system without leaving its native phase space (Hilbert space) due to peculiar 
stochastic extension of Liouville (von Neumann) equation for its statistical operator. Can one re- 
formulate the theory in terms of stochastic "Langevin equations" for its variables? Here it is shown 
that in case of classical Hamiltonian underlying dynamics the answer is principally positive, and 
general explicit method of constructing such equations is described. 
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I. Introduction. Any Langevin equations involve 
irreversibility (friction) and indeterminism (noise), as 
the classical equations which imitate interaction between 
"Brownian particle" and a fluid (see e.g. and refer- 
ences therein). Both the friction and noise represent the 
same reversible and deterministic microscopic dynam- 
ics, but usually are presumed unambiguously (additively) 
distinguishable. In general, of course, such assumption 
is wrong, because the friction itself can essentially fluc- 
tuate, as in the case of interaction between macroscopic 
vibrations of a quartz crystal and its own phonon gas (see 
e.g. and references therein). Therefore the question 
arises: how one should construct "Langevin equations" 
(interpreted loosely as a model replacement of underly- 
ing microscopic dynamics) to be sure they result quite 
accurate and thus free of artifacts? 

The answer can be formulated in the framework of "the 
stochastic representation of deterministic interactions" 
[H0000 0] at least in two widespread situations: 

i) when the dynamics is Hamiltonian while interaction 
between a system of interest, "D" , and other world, "B" , 
is described by a bilinear contribution to Hamiltonian of 
"D+B" iSli : 

H = Hd + Hb + H in t , Hi n t = J2 n D n B n ; (1) 

the marks "d" and " b" and the operators (or phase func- 
tions, in classical mechanics) D n and B n relate to "D" 
and "B", respectively; 

ii) when joint evolution operator of "D+B" , L . has 
similar bilinear form 

L = Ld + Lb + Lint , Lint = Yl n ^-n^n (2) 

The evolution operator is understood as those govern- 
ing join statistical operator, p , of "D+B" : 



dp/dt = Lp 



(3) 



In Hamiltonian dynamics, L — C(H) , where C(H) is 
quantum or classical Liouville operator, 

C(H)p=%[p,H] or C(H)p=(j%%-%%)p, 



with {q,p} being canonic variables. In case (JIJ L = 
C(H) always has the bilinear form J5J) 0,0 (the case © 
covers also non-canonic treatments of Hamiltonian dy- 
namics and, besides, essentially non-Hamiltonian and 
irreversible dynamics, and even Markovian probabilistic 
evolutions) . 

For simplicity, in this paper discussion of the Langevin 
equations will be confined by classical mechanics, more- 
over, starting from Sec. Ill, by the case only. 

II. Characteristic functionals. The statistical op- 
erator p from Eq|3 (density matrix, probability measure, 
etc.) says about current state T = Td © Tf, of "D+B" 
only. Who is interested also in its correlations with its 
prehistory, may consider one or another characteristic 
functional (CF) 



Tr d Tr b p(t, T) (exp [J £. ^(f')Q;(r(/)) dt'] 



(4) 



where Qj(T) are some phase functions (i.e. functions of 
instant system's state) and Vj(t) conjugated arbitrary 
test functions (probe functions); Tr^ and Tr^ denote 
"traces" over phase spaces of "B" and "D" , that is inte- 
grations over Tb or ; (■■■) r is conditional statistical 
average under given present state T — T(t) , and the 
right-hand side retells the left from viewpoint of exterior 
observers. Particularly, in case of deterministic dynam- 
ics the conditional averaging degenerates into replacing 
r(i') by strictly definite function of V = T(t) . 

In any case, if readdressing symbol p to the whole 
expression under the traces in Q, one can write 
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,ex P [J t>t , Y,j VjWQiP) dt'})= Tr, Tr b p , (5) 

where now, obviously, p obeys the equation 

dp/dt = {J2 j v j (t)Q j (T)+L}p (6) 

instead of Thus one reduces CF to slightly modified 
evolution equation. In fact that is a sort of famous re- 
lations between path integrals and differential equations, 
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like the Feynman-Kac formulas jToL ITU . Nevertheless, 
we once more accented the transition from (QJ to ©-© 
(see also Sec. 2 in @) because, curiously, some referees 
are not familiar with such possibility (by the way, some 
similar old examples can be found in [12ul3| ). 

III. Stochastic representation. Consider partial 
probability measure of "D" 's states, pa = Tr b p , where o 
satisfies the evolution equation ©. According to [J, [a, Hi 
if once p was factored, then later pd can be represented 
as the average 



Pd = ((aj» 



(7) 



of stochastic probability measure pd which obeys the 
time-local differential equation 



dpd 
dt 



Pd (8) 



with x n (t) and y n (t) being definite stochastic processes 
and ((...)) = ((...) y ) x = ((—) x )y statistical average with 
respect to them. Similarly, if the phase functions Qj 
wholly belong to "D" then their CF (0J can be repre- 
sented, in place of J5| and ©, as 

ex P [J t>t , Ej ViWQiV) dt']) = ((Ti-d Pd » , (9) 
where now pd is a solution of the stochastic equation 

^ = [vj (t)Qj + Vn (t)D n + C(H d + x n (t)D n )} p d (10) 

with the same random sources x n (t) and y n {t) (closely 
repeated indices imply summation). 

Notice that © and (JTDJ again exploit the Feynman- 
Kac type relations, now for stochastic evolution operator 
[Un(t)D n + C{Hd + x n (t)D n ) ] in place of L , and that 
next such instants will not commented. 

It is easy to see that x n (t) surrogate Hamiltonian 
perturbation, H d -> H d + x n (t)D n , of "D" by "B". 
What is for y n (t) , they enter ((SJ) and (fTTTf) like test 
functions conjugated with variables D n . Therefore one 
can say that y n (t) describe observation of "D" by "B" . 
But any thing under observation affects the observer. 
Hence, in other words, y n {t) represent an opposite ac- 
tion of "D" onto "B" . Importantly, this passes without 
self-action of "D" , which is the reason for peculiarity of 
random processes y n (t) ■ they are null by themselves 
(({yniih)— Vn k {tk))) = 0) althoiMi possess non-zero 
cross- correlations with x n (t) Hi, SB IE 13- 

Such corre- 
lations are responsible for energy dissipation in "D" and 
similar statistical effects. 

Quantitatively, full statistics of x n {t) and y n (t) is 
determined by separate evolution of "B" under pertur- 
bations of its Hamiltonian, Hb — > Hb + fntt)Bn , by 
arbitrary time- varying forces f n {t) 0, IE 0, S 111 ba 
this section, let T = Tb , and B n (t' 7 f,T) be values of 
the phase functions B n considered at time t' as func- 
tionals of the forces and functions of current "B" 's state 



r = T(t) at time t . Then characteristic functional of 
x n (t) and y n (t) is 

S{«(r), f{r)} = 



ee ((exp j t>tl [un{t')x n {t') + f n (t')y n (t')} dt')) = 

= Tr 6 p b {t, /, T) exp {f t>t , u n {t')B n {t', /, T)dt'} (11) 

with pb{t, /, r) being current "B" 's distribution function. 
Since B n (t' = t, f, T) = B n (T) , the expression under the 
trace in (|11|) . 

p b ee p b (t, f, T) exp {J t>t , u n (t')B n (t\ /, T) dt'} , 

satisfies the differential equation 

dpb/dt = [u n {t)B n +£{H b + f n (t)B n )] p b , (12) 

quite similar to I|1U|) . and CF (|11|) can be evaluated by 
solving this equation: 

5{u(r), /(r)} = Tr b p b (13) 

Variational differentiations of (|11|) produce the identities 



n 



Sf(r„ 



n.i x ( t j)U. m y( T -n 



Ti bPbitJ^UBfaJ^) 



(14) 



f=o 



clearly explaining the peculiarity of y n (t) ■ Besides, <|14[) 
shows the nullity of any cross-correlations between y n (r) 
and earlier x n (t' < t) 0, H, IS IS El > which is conse- 
quence of the causality principle (none perturbation of 
"D" by "B" can depend on future perturbations of "B" 
by "D"). 

IV. Fluctuation-dissipation relations. The phase 
volume conservation and generic time-reversal and time- 
translation symmetries of Hamiltonian mechanics re- 
sult in the Onsager reciprocity relations, Kubo for- 
mulas, fluctuation-dissipation theorems [jl and other 
"fluctuation-dissipation relations" (FDR) [H ESQ El- 

In |S general quantum FDR were reconsidered in 
terms of the stochastic representation. To exploit their 
classical limit, let us assume, without loss of general- 
ity, that (i) /n(-oo) = /„(+oo) = 0, (ii) B n are 
chosen so that their unperturbed mean values are ze- 
roes (i.e. ((x n (t))) = 0), and (iii) B n possess defi- 
nite time-reversal parities: B n (q, ~p) = e n B n (q,p) with 
e„ = ±1 . Besides, assume, with a loss of generality, 
that (iv) the past initial distribution function of "B" 
(before switching-on the "D"-"B" interaction) was the 
canonical one, oc exp(— Hb/T) , and (v) Hb is even: 
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Hb(q, —p) = Hi,(q,p) . Then the classical generating 

fdr [Urnm y ield 



~\u(t) 



/(r) =H{et*(-r), e /(-r)} 



1 df(r) 
T dr 

The same can be expressed Q by the equalities 

e n x n (-T) x x n (r) , 
tnVni-T) x y n (r) + T^ 1 dx n (r)/dT , 



(15) 
(16) 



where symbol >; means statistical equivalence. 

For example, averaging the product of two lines of ljl6|l 
taken with different arguments, it is easy to obtain such 
second-order relation: 



K*l{ T ) = ?Q4- K xx 



(t) 



T dr 

where 9(t) is the Heavyside step function, 
K^(t) = {{ Xl (r)x m (0))) , = ((^(r)j, m (0)» 



(17) 



and the causality principle is accounted for as prescribed 
by HH> . 

V. Distribution function. Come back to "D" as 
described by the Eg s 1715191 and ITUl using T = T d = {q,p} 
as notation for complete set of "D" 's variables. 

Equation JHJ can be viewed as generating equation 
for CF of variables D n in the system with Hamiltonian 
Hd + x n (t)D n . At that, as we already mentioned, y n (t) 
play the role of test functions conjugated with D n , while 
x n {t) are external forces. This picture is described by 
Hamilton equations and Liouville equation as follow: 

dT(t)/dt = -[C(H d + x n (t)D n )T](t) , (18) 
dp/dt = £(H d + x n (t)D n )p (19) 

Below, let T(t) — L(i, x, T, 9) denote solution of 
EalT^I with initial condition T(t = 9) = T . Besides, 
define p{t,x,T) be solution of Eg 1191 under condition 
p{tQ,x,T) — PdoiY) , where pdo is T's distribution at 
past initial time moment, to . Formally, to is the time 
when the "D"-"B" interaction was switched-on. Direct 
solution of lfH?|l reads 

p(t, x,r) = J 5{r - r(t, x, r , t )} p d0 (r )dr = 
= J s{r - T(t , x, r, t)} p d0 (r )dr = (20) 
= pdo(r(t ,x,r,t)) , 

where ${•■■} means delta- function in the phase space and 
T the initial state. At that, the group property of T's 
transformations from one time point to another: 

T{t', x, T(t, x, r , t ), t) = Tit 1 , x, r , t ) , (21) 



and the Liouville theorem about phase volume conserva- 
tion were taken into account. 

In these designations, solution of Eq|S] looks as 

p d - p(t,x,T) exp{jly n (t')D n (t',x,r,t)dt'} (22) 
with D n (t', x, T, t) = D n (T(t', x, T, t)) . 

Since y n (t) are null by themselves and null in con- 
junction with any earlier x n (t' < t) , while p(t,x,T) 
depends on x n (t' < t) only, and T(t',x, T,t) depend on 
x n (mm(t,t') < t" < ma,x(t,t')) only, one can replace the 
upper integration limit in 1221) by any value > t , in par- 
ticular, by oo . Then the exponent in 11' I'D transforms 
into the statistically equivalent functional 



St{x,y,T} = exp 



y n (t')D n (T(t',x,r,t))dt' 



(23) 



After this replacement, substitution of (1221 to (J7J, with 
use of identities iffi!)! and (|?TJl . yields 

p d (t,T) = ((p(t,x,T)S t {x,y,T})) = (24) 



= /(( <H r - r (*> x , r o, ^o)} Staix, y, Tq} )) pdo{To)dYo 

Alternatively, by averaging directly formal operator so- 
lution of EqlHl one obtains 



Pd 



eexp{fl£(H d )dt>} E{C(D) 7 D} 



PdO 



(25) 



where O symbolizes chronological ordering of the follow- 
ing operator expression (that is ordering with respect to 
imaginary time argument of H d and D n ) . 

VI. Fluctuation statistics. Similarly to preceding 
section, consider Ea llQI as generating equation for joint 
CF of variables Qj and D n . 

Now, express solutions of (|18|l and (|19|l through "D" 's 
state at arbitrary fixed time moment 9 which is different 
from t , that is solve l|18|) and (|19|l under initial condi- 
tion T{t' =9) =T (thus T{t' = 9,x,T,9) — T). Then 
solution of EqE3 can be implicitly formulated as 

p d (t,T(t,x,T,9)) = p(9,x,T)x 
exp { / [ Vj tf)Qi (f , x, T, 9) + y n (t')D n (t', x, T, 9)} dt'} 

(26) 

with Qj(t',x,T,e) =Qj(T(t',x,T,9)). Substituting 
to (jHJl and taking into account the Liouville theorem (the 
phase volume conservation under arbitrary Hamiltonian 
evolution), at t — > oo we have 



(exp { / Vj (i')Qj (*') dt'}) = ((Tr d exp { / v 3 {t')Qj (f , x, T, 9) dt'} p{9, x, T) Sg{x, y, r}» (27) 
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In terms of various statistical moments of variables Qj (omitting their indices) 

( Q(h)... Q(t k ) } = Tr d ({Q(h,x, T, 8)... Q(t k ,x, T, 8) p{9, x, T) S e {x, y, T} }) (28) 

In particular, if 8 — > to then p(9,x,T) turns into the initial distribution, p d o{F) , definitively independent on x n (t) : 

(Q(t 1 )...Q(t k )} =Tr d p do (T){(Q(t 1 ,x,T,to)...Q(t k ,x,T,to)S to {x,y,T})) (29) 

with r representing the initial state To ■ The same expression results from (12811 after substitution of (121 H and j|21(l . 
Alternatively, quite similarly to l|25|) . 

(exp {/ ViWQtf) dt'}) = Tr d § exp{ + C(H d )] dt'} 3{C(D), D} Pd0 (30) 

The functional 5 here, defined by at once accumulates all information about "B" which must be used when 
evaluating fZQ-fgft. 



VII. Self-interaction through environment and 
"scattering operator". It is useful to empha- 
size rather interesting resemblance between Eo.12%1 
or Ea l29l and expressions for scattering amplitudes, 
Green functions, etc., in quantum theory of fields 
and many-particle systems (see e.g. 0, Il7|'l. 
If draw an analogy from Qj and x n (t) to elec- 
tron operators and radiation field, respectively, then 
the averages ({Q(ti,x,T,to)...Q(t k ,x, T, to))) corre- 
spond to lowest-order perturbation approximation, while 
{(Q(t 1 ,x,r,to)...Q(t k ,x,r,to)S to {x,y,r})) in EqEl 
exactly summarizes all the orders of "D"'s interaction 
with its environment. The analogy continues in that the 
"complete multiple scattering operator" St {x,y,T} by 
itself behaves like unity: 

({s t0 {x >y ,r})) =i (31) 

This identity clearly follows from EqEl at t — » to 
and is easy explainable if notice that in any term of 
St {x, y, r}'s series expansion over y n (t') and x n {t') 
most late time argument belongs to some of y's. 

According to (fTT|l and (fT4"| , separately x n (t) are noth- 
ing but noise of free unperturbed environment, like "zero, 
or vacuum, fluctuations". However, along with y n {t) in 



St {x,y,T} they represent actual noise of the environ- 
ment, including its directional response to the system's 
motion, in the form of both renormalization of primor- 
dial "D" 's dynamical properties and appearance of new 
ones: relaxation, "spectral lines broadening" , etc. 

VIII. State-dependent noise and the fiction of 
friction. In jj.Lj.Lj the words "Langevin equation" were 
addressed to objects like JHJ) or (fTT)f) which emerged as 
stochastic extensions of the Liouville equation for proba- 
bility measure of T . In usual sense, Langevin equations 
must be a stochastic extension of the Hamilton equations 
for r themselves. Besides, one would want these equa- 
tions to involve some "realistic" noises only but not aux- 
iliary "ghost" noises like y n (t) . The latter requirement 
means that desirable equations are certainly not literal 
consequence of the basic Eqs|S|and[TO] Instead, Langevin 
equations must be especially constructed as their exact 
statistical equivalent (or at least close approximate one). 

In should be underlined that, at such target setting, a 
"size" of system "B" is insignificant (no matter e.g. is a 
Brownian particle macroscopic or as small as molecules) . 

With the formulated purpose, let us return to Eal2Hl 
choosing arbitrary functions Qj(T) as delta- functions 
S{T — 7} and their index as time. Then Eg 1291 produces 



^{7} = ( IL S{T(t) - 7 (t)} ) = / (( n t 6{T(t, x, To, to) - 7 (t)} exp {/ y(t)D( 7 (t)) dt} )) Pd o(To)dTo (g2) 

which represents probability density functional for the whole system's trajectory. Here all non-principal indices are 
omitted, and the delta-functions have allowed to replace D(T(t, x, To, to)) in the exponent by D(^(t)) . 

The simplest construction of Langevin equations follows directly from careful "visual" investigation of Ea l32l This 
shows that Eq|221can be rewritten as 

^{7} = / «n* S{T(t, z, T , t ) - 7 (i)} }) 7 Pdo(r )dT Q , (33) 

with new random forces z n (t) in place of x n (t) , if conditional statistics of z n (t) is defined by formulas 

((z(t 1 )...z(t k ))y ((x(tx)...xQk k ) e^ify^D^dt) )) (34) 
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The brackets (( ... )) 7 here have the sense of conditional averaging under given system's trajectory j(t) . 

At that, the role of Langevin equations governing the variables T(t) = T(t, z,T, to) and Q(t) = Q(T(t, z,T,to)) 
belongs to nothing but merely the Hamilton equations: 

dT(t)/dt = -[C(H d + z n (t)D n )T](t) (35) 



Notice that in view of identity (|31jl the averaging pro- 
cedure defined in l|34fl automatically satisfies the nor- 
malization condition ((l)} 7 = 1. Besides, due to the 
above mentioned statistical peculiarities of y(i)'s, result 
of the averaging always agrees with the causality princi- 
ple: the moments <|34|) in fact can depend on ^(t) with 
t < max(ii, ...,t k ) only. 

Formally, the two above expressions, (|3~4fl and (J35), 
already define what can be named "exactly equivalent 
Langevinian form of the stochastic representation" . It 
clearly emphasizes statistical nature of dissipation and 
friction: even if being present they still are hidden in- 
side z(t)'s statistics, as cross-correlations of x(i)'s with 
y(i)'s. To see them evidently in (|35[) . we have to with- 
draw them from 134(1 in some reasonable approximation, 



with corresponding redefinition of noises z n (t) . 

IX. Unbiased noise and Langevin equations. 

With the above pointed purpose, first, assume, naturally 
and without loss of generality, that (x) = . Then de- 
sired dissipative contributions to (|35|l . together with the 
renormalization corrections of non-dissipative terms, can 
be identified among mean values of z(t)'s. 

Second, consider cumulants (semi-invariants) n a f3 = 
((xi, ...,x a , yi, yp )) . For brevity, here the subscripts 
unify indices and time, and commas do emphasize that 
comma-separated multipliers arc subject to purely irre- 
ducible correlation of (a + /3)-th order ("Malakhov's cu- 
mulant brackets"). Then CF can be symbolically 
written as 



= exp[K{u, /}] = exp 



a=2 a, (3=1 ' 



(according to our assumption, k\o = (x) =0). Decompose it into two multipliers: 



s {w, /} = -{u, /} /} , ^{w,/} = exp 



3=1 



, E{u, /} = exp 



OO OO 



X! X! Ka i 

a=2 0=0 



' a!/?! 



(36) 



Correspondingly to this factorization of CF, both the original noises, x(t)'s and j/(t)'s, divides into two components: 
x = x + x and y = y + y , where two pairs {x, y} and {x, y} are mutually statistically independent. 

It is easy to prove that for arbitrary functional <E>(x) and arbitrary function / the equality holds as follows: 



«$(x)exp(/y)))=$(X(/)) , X(f) = Y,"ipf/ffl 

13=1 



(37) 



This is because the pair {x, y} , in accordance with 114JI . describes merely conditional mean value of "B"'s response 
to its perturbation by forces / . Applying the decomposition (|36|) to l|32l) , with the help of (|37H we obtain 



W{l} = f((Ut 6 i r ^ x (D(l)) + x, T , to) - 7 (t)} exp {/ y(t)D(>y(t)) dt) )) p d0 (T ) dT Q 
The mean response X(f) , defined by l|37|) . with / = D(p() , after restoration of its temporal index, reads 



(38) 



X(t,f) = ^ln~{uJ}\ u=0 = J \{x{t) 1 y{t 1 )))f{t 1 )dt 1 + l - J J ((x(t), y(h), y{t 2 ))) f {h) f (t 2 )dhdt 2 + ... , (39) 

where because of (|14|) all integrals are in fact taken over tj < t . It is useful to notice also that, due to the causality, 
Jacobian of mutual transformations between 7 and T is unit. 

Scanning (|38|l in comparison with ll-il't and l|33|l , one evidently comes to another form of the probability functional: 
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W{j} = / ((lit H r(t, X(D(j)) + z,r , h) - j(t)} )) 7 p d0 (T )dT 



(40) 



where statistics of renormalized (in fact merely biased) noises z(t) is now described by 



«*(*!)...*(**) » 7 = ({x(t 1 )...x(t k )exp{fy(t)D(j(t))dt})) 



(41) 



At that, correspondingly to (|40|) . r(i) = T(t, X(D(T)) + 2, To, to) , that is stochastic Hamilton equations l|35|) change 
to the stochastic integro-differential equations 



dT(t)/dt = - {jC(H d )T](t) - (X n (t,D(r)) +I n (t))[C(D n )T}(t) 

I 



(42) 



As prescribed by (|36|l and (|41() . here the noises z n (t) 
have certainly zero mean values, while any dissipa- 
tive effects of interaction with "B" are separated in 
X n (t, D(T)) . Hence, Eqs02l can be by now enough surely 
named "Langevin equations" . 

X. Discussion. Of course, the above result is rather 
trivial one. However, from the point of view of appli- 
cations and practical computability, it is not quite sat- 
isfactory. The matter is that numeric modeling of noise 
essentially conditioned by the system it drives is generally 
difficult task. It would be better if the noise was reduced 
to unconditioned random quantities, for example, 



z(t) = z W(t) + f zW (f,ti)/(ti) dti ■ 
+ i fj z^\t,hM) f(ti)f(t 2 )dtidt 2 + 



(43) 



where z^°'(t) = x(t) = x{t) is unperturbed noise, 
z^(t,ti) represents stochastic linear response of "B" to 
its perturbation, etc., and all z^ are some zero-average 
random functions independent on the forces. In par- 
ticular, z^(t,ti) includes fluctuations in linear friction 
(whose average was contained in first term of (|2HJ))- 

The only situation when Eqs021 finalize the analysis 
is when the noises z(t) are state- independent, that is 
z ( n ) = o for all n > . But this is unlikely realistic 
situation since in general it is forbidden by restrictions 
which follow from the phase volume conservation and mi- 
croscopic reversibility. For concreteness, if "B" is equilib- 
rium thermal bath (thermost at), these restrictions are ex- 
pressed by FDR <H3 [II El GJ, [13 or equivalently JTSJ 
(notice that FDR for internally non-equilibrium baths 
also were considered in 0, 0). If noises z(t) are in- 
deed state- independent, this means that n a p — for all 
a > 2 and (3 > 1 . Then the second row from (|16|1 clearly 
implies that in such case the equalities k q o = also 
should hold for all a > 3 . In other words, the noise z(t) 
can be purely state-independent only when it is purely 
Gaussian. Moreover, then the same FDR prescribe that 
= for all (3 > 2 , that is average response of "B" is 
purely linear. 

Thus we come to the trite "linear Gaussian thermo- 



stat" when Easl41l and 1391 reduce to 

X n (t,D(T)) = -K*? n (0)D m (T(t))/T + 



(44) 



T 



K™(t-t>)—D m (T(t>))dt' 



((z n (h)Zm(t2) = K™(ti-t 2 ) 

Here FDR (|17|) is used, and it is taken in mind that all 
higher-order cumulants of z(t) are zeros. Discussion of 
more interest models will be done elsewhere. 

XI. Example: oscillator. Consider nonlinear oscil- 
lator, assuming that "B" is "linear Gaussian thermostat" 
while interaction with it realizes in potential way through 
two statistically independent channels as follow: 

H d = p 2 /2m + U (q) , D^T) = -q , 
D 2 (T) = q 2 /2 , K& = K%=0 

The first channel corresponds to usual thermal excita- 
tion, and the second to thermal parametric fluctuations 
in frequency of oscillations. The Eas 1421 and 1441 yield 

dq{t)/dt = p(t)/m , 

dp(t)/dt = -dU(q(t))/dq(t) +5i(i) +x 2 {t)q(t)- 
- £«, *>(*WT - 

- q(t) fL K 22 x (t - 3(*>(f) dt'/T , 

where v(t) = dq(t)/dt is velocity, x n {t) are mutually 
independent normal random processes, K** are their 
correlators, and 

U(q) = U (q) - ifff(0) q 2 /2T - K% (0)q 4 /8T 

is renormalized potential. Hence, correspondingly, there 
are two channels of friction and dissipation, and the fric- 
tion channel conjugated with thermal parametric fluctu- 
ations is essentially nonlinear. Similar examples concern- 
ing thermal fluctuations in capacities of electric circuits 
were considered in [Tsl ]. 
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XII. Conclusion. For particular variant of the 
"stochastic representation of deterministic interactions" 
concerning classical Hamiltonian mechanics, we have 
demonstrated that by request it can be completely re- 
formulated in terms of "Langevin equations" for internal 
variables of an open system. These equations are wholly 
housed in its own phase space and are free of the peculiar 
auxiliary noises y n (t) , distinctive for initial "stochastic 
representation". At the same time, y n (t) remain useful 
undercover instrument, being responsible for conditional 
statistical dependence of actual noise on trajectory of the 
system driven by it. 

This Langevinian form of the theory seems more vivid, 
although, probably, it will occur less appropriate for prac- 
tical analysis of complicated noise statistics. Besides, the 



original "Liouvillian form" at once covers quantum me- 
chanics as well. 

What is for its quantum Langevinian equivalent, still 
it remains unexplored. Notice that quantum Langevin 
equations for important special case of Gaussian lin- 
ear thermostat were exhaustively considered in 0]. Of 
course, more general situations also were under many 
considerations (see e.g. pofl. 

But recall that the question under our principal and 
pragmatic interest is how much non-Gaussian non-linear 
generalization of quantum Langevin equations can be de- 
veloped if do it wholly within native Hilbert space of an 
open system under consideration and with use of com- 
mutative ( c-number valued) noise sources only. 
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